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Abstract 

The grand canonical partition functions for primitive droplet nucleation mod- 
els with an excess energy e„ = —fin + crn^"'', r] = 1/d, for droplets of n con- 
stituents in d dimensions are calculated exactly in closed form in the cases 
d = 2 and 3 for all (complex) fi by exploiting the fact that the partition 
functions obey simple linear PDE. 
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The picture that Ist-order phase transitions are initiated by the formation 
(homogeneous nucleation) of expanding droplets of the new phase within the 
old phase is a familiar and important one (see the reviews [1-13]). In its most 
primitive form the droplets are considered to consist of n constituents (e.g. 
droplets of Ising spins on a lattice or liquid droplets of molecules etc.), the 
"excess" energy of which is given by a "bulk" term proportional to n and 
a "surface" term proportional to n^~^,ri = 1/d, where d > 2 is the dimension 
of the system: 

Cn = —p^n + an^~^ , i] = 1/d . (1) 

In the case of negative Ising spin droplets formed in a background of positive 
spins by turning an external magnetic field H slowly negative, below the 
critical temperature, the coefficient fi in Eq. (1) takes the form fi = —2H. 
For liquid droplets of n molecules condensing from a supersaturated vapour 
one has ft = fi — fic, where is the chemical potential and Hc its critical value 
at condensation point. 

Assuming the average number u{n) of droplets with n elements to be given 
by a Boltzmann factor, 

uin) (X e-(^'- ^ f3 = -L. ^ (2) 

Kb! 

and that the droplets form a noninteracting dilute gas leads to the grand 
canonical potential ipd per spin or per volume 

n=0 

t = Pfi, X = Pa ; p : pressure , 

d^jjd = -Udp + ndt . (4) 

(Physical reasons may require to start the sum (3) not at n = but at some 
finite rio > 0. This does not affect the following conclusions and can easily 
be taken care of. It is mathematically convenient to start at n = 0.) 
The series (3) converges for t < only, which follows e.g. from the Maclaurin- 



Cauchy integral criterium [|14|. In applications to metastable systems one is 
interested, however, in the behaviour of ipd{t,x) for t > which calls for an 
analytic continuation in t or in the fugacity z = [|15], 0, |l^ . 



The series (3) has recently been discussed in connection with the canonical 
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quantum statistics of Schwarzschild black holes |T^, |T^, |20 
Notice that it obeys the linear PDE 

df-'^, = i-iy ^t^|J, , (5) 

in particular 

dttp2 = dlij2 for ci = 2 , (6) 

d^ip3 = -dliJs ford = 3, (7) 

which will be helpful to find ipd for t > 0, in the following especially for d = 3. 
Before discussing the exact continuation of 1^2 and ip3 into the complex t- or 
2-plane, let me recall the known saddle point approximations for ijjd if n is 



turned into a continuous variable and t > [16 



V;,= / cine^^/^'^-^^'"') . (8) 



Setting n = {a / fi) u yields 

^ ) r rfMM^"ie^^(") , hiu) = a(u'' - u''-^) , a = . (9) 
fi J Jo /i" ^ 

The function h{u) has extrema at m = = and ai u = ui = 1 — rj 
with h{ui) = —ar]{l — r]Y~^, h"{ui) = a{d — 1)(1 — . According to 

the standard analysis the integral (9) has for large /3 the asymptotic 
expansion 

d/2 



i>d-il- riY^sj^^ (^j e-/^«^(l - ^)'"' {^ + 0(l//3)) . (10) 

Here the path in the complex w-plane goes from mq = to Ui and then parallel 
to the imaginary axis to +ioo. Thus, only half of the associated Gaussian 
integral along the steepest descents contributes! The leading terms in (10) 
are purely imaginary! For d = 2 and 3 (here see also Ref. [16]) one has 

4-\t,x) = ^^e-V(4t), (n) 



.2v/?x3/2 



4 X' 
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^t\t,x) = ^-^^e 27 . (12) 
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Whereas ijj2°°^ is a solution of the heat equation (6), tjj^^ is not an exact 
solution of the corresponding equation (7)! We shall see below how this is to 
be understood in terms of the exact ip3{t, x). 

For d = 2 the function il)2{ti has been determined exactly for t > in closed 
form in Ref. |T^ by following Lerch's observation that the relation 



a/tt io 



^' e-^V(4r) - nr ^. r ^^^^^^ ^y-nr ^^3) 



2^h T^l^ 

converts the series(3) into a geometrical one under the integral sign. Here 
k{t, x) = — ^e-^V (4^) = -2d.,K{t, x) , (14) 



where 



i^(t,a;) = ^e-^V(4t) (15) 



is the "heat kernel" which obeys the heat equation (6). So does K{t,x). 
Now the series (3) can be summed exactly for t < (x > 0) and then 
continued: 

Mz = e\x) = / dTK{T,x) J- r (16) 

Jo i_e[t-T) 

X roo g-xV(41nM'l ^ 



^ ^ I du ^7:5 . (17) 

2^ Ji In^l^u u-z ^ ^ 

The last relation shows that ip2{z,x) can be continued analytically into the 
whole complex 2;-plane except for a cut from 2; = 1 to 2; = 00. 
The discontinuity of 1^2 across the cut {z real and > 1) is given by 

lim [7/^2(2 + ie) - ^2{z - it)] = 2mk{t, x) , (18) 
e ^ 0+ 

and, if one approaches the cut from above, the limit 

lim 'ip2{z + ie^x), 2; real and >1, 
e ^ 0+ 
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is no longer a real-valued function of but has a nonvanisliing imaginary 
part 

^Mt, x)] = nK{t, x) = ^e-^V(4t) , (19) 

which agrees exactly with that of the saddle point approximation (11)! 
The real part ^[ip2{t,x)] is given by the principal value integral 

nMt.x)] = p.v. / dTKiT,x)- — . (20) 

Jo I — ^ 

In order to prepare for the method by which to determine i/'s, let us arrive 
at the above result for il)2 in a different way: 

As the heat equation (6) is invariant under the scale transformation t — »• 
At,x A^x, A > 0, it has solutions of the type giy = x'^/t), where g{y) 
obeys the ODE 

9"+i\ + Yy^9' = 0. (21) 



The solutions (up to additive and multiplicative constants) are (see Ref. |23[| ) 

g{y) = y-'/'e-y/^wik = -1/4, m = 1/4; y/A) , (22) 

where w is any solution of Whittaker's standard form of the confluent hy- 
pergeometric equation |2^. As we want to associate the solution g with the 
imaginary part '^[ip2] in (11), we take w{— 1/4, 1/4; y/A) = aW_i i{y/4),a: 
constant, where Wkm is Whittaker's "W"-function which has the property 
[11 

, 1 /o 2 /"CO 2 

u-y^e-^/^^ W_i^i{u^) = 2j dve-"" =2Eifc{u). (23) 
Here Erfc(u) is the complementary error function. Thus we get 

g{y) = const. Eifci^^) . (24) 

As 

a,Erfc(i V^) = ^e-^V(4t) = ^K{t, x) , (25) 

we see that we may arrive at the solution (16) by essentially folding the 
imaginary part dtg{y) with the function 1/(1 — e*^"^). (Note: If F{t,x) is a 
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solution of the heat equation, then any derivative or integral of F{t, x) with 
respect to t or x is a solution, too!) 

The normalization of g{y) is fixed by observing that for A > we have (by 
rescaling r A^r) 



^2(t,Ax)=/ dTK{T,x)- (26) 

Jo 1 — exp(t — At) 

and that therefore 

lim ^'2(^5 Ax) = T , lim ■ilj2it,Xx) = 1 . (27) 

A— »0 X — gl- A— >oo 

We are now ready to proceed in the same way for ip^: As the PDE (7) is 
invariant under the scaling [|19[ t ^ \t,x A^^^x, the ansatz g{y = x^/t^) 
leads to the ODE 

^'+4 + V + |(- + Ay = 0, (28) 
21 y 9 y y' 



which has solutions (see Ref. p3| ) 



2 

^'(y) = r'e~27\(l/2,l/6;^y) . (29) 



The choice w = Wi i{^y) provides the desired result: 
With 

2 



we have 



2 

%(l/) = ^e"27V.,.(^l/), (31) 

which for large i.e. large x or small t, takes the asymptotic form 

4 

4 x2/3 xV^^ 

27 . (32) 
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Except for a factor v^37r/2 this is the same as i)^^^ in Eq. (12). 
Let us, therefore, try the ansatz 

ip3{t,x) = / dTK:i{T,x)- ^— - , (33) 

Jo 1 — e*- ^ 

2 3/2 

r 2 '5 ^ (' 

where a is a normahzation constant to be determined. 
That the function ^/'s of Eq. (33) is indeed the right one can be seen as follows: 
It is a solution of Eq. (7), because by construction i^3(t, x) is a solution of 
that equation (the partial derivative of ■j/'s with respect to t can be replaced 
under the integral by the negative derivative with respect to r, followed by 
a partial integration with repect to r and observing that K^^ vanishes at the 
boundaries, see Eqs. (32) and (39)). 

Furthermore, the imaginary part of coincides with that of Eq. (12). Fi- 
nally, the required boundary conditions can be fulfilled by an appropriate 
choice of a: 

Rescaling r X^/'^t yields 

7/>3(t,Ax) = / drK^ir^x)- . (34) 

In order to have (see Eq. (3)) 

lini?/^3(t, Xx) = ^ . , lim ip^it, Ax) = 1 , (35) 

A^O I — g'' A^oo 

we need 

/•LXJ 

dTK3{T,x) = 1 . (36) 



Because the relation 

r%-^W..i2u)= ^,, =2,/^ (37) 
Jo u 2'6^ ^ cos(7r/6) V3 ^ ^ 



holds ||25|, the normalization (36) requires a = \j'i/Ti, so that finally 

9 

— 3/2 
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Let me add that Ks{t,x) has no zeros on the positive real r-axis for x > 
[24] and therefore is strictly positive there. Furthermore, for u ^ one has 

m 



W.,.{v) ^ ^nV3(i + 0{u)) + ^^/^'il + 0{u)) , (39) 
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which is of interest for the behaviour of Ksit, x) in the limits x — or t— >■ oo. 
One can check the result (38) by expanding the factor 1/(1 — e'-*"'^^) in the 
integrand of (33) for t < in a geometrical series and comparing the coeffi- 
cient 

/(x, n) = \^ r dr-e-^^"'^" Wi i {^x^'/r^) e"^ ^ (40) 

V TT JO T 2'6 27 

of exp(nt) with the one - exp(— xn^/^) - of the original series (3): 
Because /(Ax, n) = /(x, A^/^n), we have /(x, n) = g{xn'^^^) and it is suffi- 
cient to consider (40) for x = 1. f{n) = /(x = l,n) can be interpreted as 
the Laplace transform of 

1 2 4 

h{r') = TT^e-iTT. ' = I dre-^^h{r') ■ (41) 



In that case one has (see Ref. K 



/H = ^r^e-^'«V4^(lM^ (42) 

VTT Jo 

where g{y) is the Laplace transform of h{t = r^). 
As[i3 

g{y) = dte-yh{t) = -i= ^^1/3(4^^) , (43) 

where Ky{z) is the modified Hankel function |28|, we finally get from Eqs. 
(42) and (43) the desired result (compare Ref. p9|): 

As (see Ref. [13) 



'2z z 
W^i,i(2^) = ^{W,^.{2z) + W,^.{2z)) = —{K._{z) + K.{z)) , (45) 
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K3{t,x) may also be expressed by these special functions. 
Having determined the exact partition functions, one can now study the asso- 
ciated thermodynamical properties, e.g. the behaviour of the magnetization 
which is essentially given by the derivative of ipa with respect to t, the de- 
tails of the metastabilities etc. However, improving the mathematics does 
not eliminate the many shortcomings of the model concerning the physics it 
is supposed to approximate (see the Refs. [1-13])! Nevertheless, it is always 
pleasing to have exact solutions for nontrivial models. 

I thank G. Rocpstorff and T. Strobl for a critical reading of the manuscript 
and the former for drawing my attention to Ref. [17]. 
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